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(54) Verf ahren zur Infflalisierung und Aktualisierung elnes Netzmodells 

durch nachstehende Schritte gebikJet: Bereitstellung 
erfaBter Netztopologiedaten in Form einer Terminal- 
Adjazenz-Matrix, Umformung der Matrix in eine obere 
Dreiecksmatrix, Initiaiisierung der oberen Dreiecksma- 
trix, Bildung der Matrixprodukte fOr Zeilen der oberen 
Dreiecksmatrix, Bildung einer Pfad-Tabelle ausder obe- 
ren Dreiecksmatrix Aktualisierung der oberen Dreiecks- 
matrix im Fall einer Anderung der erfaBten 
Netztopologiedaten, Aktualisierung der Pfad-Tabelle 
und Kennzeichnung der Komponentendes Netzgraphen 
anhand der Pfad-Tabelle. 



(57) Die Erf indung bezieht sich auf ein Verfahren zur 
Initialisierung und Aktualisierung eines Netzmodells, 
wobei die Netztopologie in Form eines Netzgraphen 
gebildet wird und die Aufgabe zugrundeliegt, ein Verfah- 
ren anzugeben, das mit einem Rechnerprogramm reali- 
sierbar ist das mit einer kurzen Rechenzeit abarbeitbar 
ist. Solche Verfahren finden beispieisweise in Netzleit- 
rechnern eines elektrischen Energieversorgungsnetzes 
Anwendung. wobei die ermittelte Netztopologie Basis fQr 
Anzeigen oder Berechnungen ist. 

Beim erfindungsgemaBen Verfahren wird die Netz- 
topologie unter Anwendung und Erweiterung von 
Sparse-Vector-Verfahren und Sparse-Matrix-Verfahren 
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Beschreibung 

Die Erfindung bezieht sich auf ein Verfahren zur Initialisierung und Aktualisierung eines Netzmodells, wobei die 
Netztopologie in Form eines Netzgraphen gebildet wird. 
s Das Netzmodell kann sich auf Netzwerke unterschiedlichster Art beziehen, z.B. auf elektrische Netze, Gasnetze 

oder Nachrichtennetze. Nachstehend wird jedoch zur Vereinfachung stets von einem elektrischen Energieversorgungs- 
netz gesprochen. 

Zur Uberwachung und FQhrung elektrischer Netze ist die Kenntnis und Darstellung des topologischen Netzzustan- 
des sowie die Verfolgung der Netzspannung von grundlegender Bedeutung. In den Leitrechnern von NetzfQhrungssy- 
w stemen werden deshalb topologische Netzmodelle als ProzeBabbilder des zu fuhrenden Obertragungs- oder 
Verteilnetzes aufgrund statischer Informationen, den Topologie-Daten des Netzes initialisiert und gehatten und aufgrund 
von ferngemeideten Schaltermeidungen und MeBwerten (dynamische Informationen) foiHaufend aktualisiert 

Die topologischen Netzmodelle werden fur unterschiedliche Zwecke verwendet, z. B. zur graphischen Darstellung 
des aktuellen Netzzustandes, vorzugsweise durch Einfarbung von Teilen des dargestellten Netzes, wie z.B. beschrieben 
is von W. Grein et al. in "Dynamic Network Colouring", Proceedings of the Eighth Power Systems Computation Conference 
(PSCC), Helsinki, 19-24 August 1984, Seiten 489 bis 496. Auch als Auskunftsfunktion fur globale Verriegelungen, z.B. 
Schalten auf Erde, oder als Grundlagefur Berechnungen. Solche typischen Berechnungen kflnnen sich auf einen kon- 
sistenten Netzzustand (State Estimation} oder simulierten Netzzustand (LastfluBberechnung) beziehen. Weiterhin gibt 
es separate Berechnungsfunktionen zur Bestimmung von lokalen Feldtopologien. 

20 Bekannte Verfahren zur Initialisierung und Aktualisierung beruhen auf Aigorithmen fur Graphen oder sparlichen 
Adjazenz-Matrizen. Diese Verfahren machen es erforderlich, um den Netzzusammenhang zu ermitteln (Initialisierung), 
den Graphen des Netzes zu durchsuchen und fur jede Komponente einen Spannbaum zu bestimmen. Fur groBe Netze 
bedeutet dies einen erheblichen Aufwand. Deshalb wird bei einer Schalteranderung versucht, nur die wirWich notwen- 
digen Anderungen (Update) an den Datenstrukturen, die den Netzgraphen oder die sparliche Adjazenzmatrix beschrei- 

25 ben, durchzufuhren. Die Bestimmung dieser notwendigen Anderungen kann je nach Netzanderung sehr aufwendig sein. 
Das Verhaitnis von Update zu Initialisierung ist fur merWiche Anderungen deshalb in der GrOBenordnung von 1 :10. Im 
Stdrungsfall laufen in das Leitsystem eine Vielzahl von Schaltermeidungen ein, die sequent ell von der Topologiefunktion 
verarbeitet werden mQssen. Um den Leitrechner in kritischen Situation nicht mit einer fortlaufenden sequentiellen Ver- 
arbeitung topologischer Updates zu Qberlasten, werden Oberlastverfahren eingesetzt, die je nach Situation einen Topo- 

30 logie-Update oder eine Topologie-lnitialisierung veranlassen. Alle bisher bekannten Update-Verfahren sind aufwendg, 
komplex und im mathematischen Sinne nicht abgeschlossen und deshalb fehleranfallig. Aufgrund der Komplexrtat sind 
die entwickelten Programmpakete sehr umfangreich sowie die damrt verbundene Programmpf lege und die Einarbeitung 
neuer Mitarbeiter in die Anwendung diese Programme sehr aufwendig. 

Auf dem Gebiet der Sparse- Vektor- Verfahren und -Matrix- Verfahren sind Aigorithmen entwickelt worden, die auch 

35 fur groBe elektrische Netze im Rahmen von Netzberechnungen einsetzbar sind. Solche Verfahren sind z. B. beschrieben 
von V. Brandwajn et al. in "Sparse Vector Methods", IEEE Transactions on Power Apparatus and Systems, Vol. PAS- 
104, No. 4, February 1985, Seiten 295 bis 301. In dieser VerOffentlichung sind auch Begriffe wie Netzgraph Oder Pfad- 
Tabelle eriautert, die in dieser Beschreibung verwendet werden. Die von Brandwajn et al. beschriebenen Berechnungs- 
verfahren setzen eine Netztopologie bereits als vorhanden voraus. Eine Untersuchung zu AnwendungsmOglichkeiten 

40 solcher Verfahren zur Topologiebestimmung ist nicht bekannt. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Bestimmung der Netztopologie anzugeben, das eine 
vergleichsweise wenig aufwendige und schnell durchfuhrbare Initialisierung und Aktualisierung eines Netzmodells 
ermOglicht. 

Diese Aufgabe wird durch ein Verfahren zur Initialisierung und Aktualisierung eines Netzmodells gelOst, das die im 
45 Anspruch 1 angegebenen Verfahrensschritte aufweist. 

Vorteilhafte Ausgestaltungen sind in weiteren Anspruchen angegeben und der unten stehenden Erlauterung der 
Erfindung anhand von Beispielen zu entnehmen. 

Das erfindungsgemaBe Verfahren ermdglicht eine um GrOBenordnungen schnellere Initialisierung oder Aktualisie- 
rung eines Topologiemodells als bekannte Verfahren. AuBerdem ist die Aktualisierungs-(Update-)Zeit weitgehend unab- 
so hangig von der NetzgrGSe. 

GemaS einer vorteilhaften Ausgestaltung wird mit einer topologischen Adresse gearbertet, die die Dateneingabe 
vereinfacht, und die es ermdglicht, lokale topologische Funktionen, z.B. Feldtopologien in die Netztopologieberechnung 
zu integrieren. MeBwerte werden dabei auch als topologische Objekte modelliert wodurch sich die Dateneingabe ver- 
einfacht und eine exakte Beschreibung auch solcher Netzkonfigurationen mdglich ist, bei denen bisher zusatzliche 
55 Pseudo-Netzelemente vorzusehen sind. 

Netzberechnungen, die auf der modellierten Netztopologie basieren kdnnen durch Anwendung des erfindungsge- 
mSBen Verfahrens wesentlich beschleunigt werden, weil notwendige Datenstrukturen und Veraibeitungsschritte schon 
in dem vorgeschlagenen Verfahren enthalten sind. Neue Impulse und Aspekte ergeben sich auch fur Programme wie 
die (n-1 )-Ausfallrechnung und die topologische Beobachtbarkeitsanalyse. Alle diese Programme haben die Auswertung 
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von ToDoloaievarianten des Netzes gemeinsam. Sie prof itieren von dem schnellen Update-Verfahren Deshalb werden 

bei denen Auswirkungen von NetzSnderungen sehr schnell festgesteHt werden mGssen, 
wLl^ erwartet. Das k6nnen z.B. Topologieverfahren fQr Nachrichtennetze Oder auch Opfm*- 

rungsverfahren mit Hilfe von Network-Flow Algorithmen sein. 
s Zur Eriauterung der Erfindung herangezogene Zeichnungsfiguren sind nachstehend aufgelistet. 

Eszeigen: 

Fig. 1 Stationen und EVU-Stationen, 
Fig. 2 Aulbau eines Feldes in einer EVU-Station. 
10 Fig. 3 Adressentransformation, 
Fig. 4 MeBwertzuordnung, 
Fig. 5 Beispielgraph, 

Fig. 6 Terminal-Adjazenz-Matrix und deren LDU-Zerlegung, 

Fig. 7 Initialisierung der qberen Dreiecksmatrix 
is Fig. 8 Ganzzahlig initialisierte obere Dreiecksmatrix. 

Fig 9 Pfad-Tabelle und Pfad-Graph (Beispiel), 

Fig. 10 Mathematische Formulierung von Initialisierung und Update, 

Fig 11 Ganzzahlig aktualisierte obere Dreiecksmatrix, 

Fig. 12 Pfad-Tabelle und Pfad-Graph fQr modifizierten Beispielgraphen, 
20 Fig. 13 Knoten-Zweig-Netz, 

Fig. 14 Ableitung der Knoten-Admittanz-Matrix, 

Fia 15 Testgraphen for Leistungsmessung, 

Fig 1 6 Fiktive Zweige in Abhangigkeit der NetzgrGBe und des Vermaschungsgrads, 
Fig 1 7 CPU-Zeit zum Aufbau der statischen Topologie (Testbeispiel), 
25 Fig. 1 8 CPU-Zeit zur Initialisierung der statischen Topologie (Testbeispiel), 
Fig. 1 9 CPU-Zeit fQr die Aktualisierung (Testbeispiel). 

Die unten stehenden Beispiele zur Eriauterung des Verfahrens basieren auf der nachstehenden vorteilhaften Def i- 
nition iT^er^ intemen topologischen Adresse fQr elektrische Obertragungs- und Verte. netze. Diese Defin.- 
tT^ S erm.glicht es, lokale topologische Funktionen (Feldtopologieberechnung) ,n d.e 

Verbindungen *ri hierbei impedanzbehaftete Netzelemente ^wie Leitung, 
S^M^ «*■- ^ anzlose Netzelemente wie Leistungsschalter, Tr^ner oderl M£r 
S^H^ sind externe topologische Adressen im Netz urxJ kOnnen Netzelementen wie Sarmiel ^»n^ 
sSen Ser Teilen von Netzelementen z.B. den beiden Enden (Terming) einer Latung zugeordnet werden. FQr 
elektrische Obertragungs- und Verteilnetze wird nun folgende topologische Hierarchie definiert. 



30 



35 



40 



45 



1 Station (oder Geographische Station) 

2 EVU-Station (oder Spannungseb^ne) 

3 Feld 

4 Terminal 



so 



55 



Fin 1 zeiqt beispielhaft ein Net. bestehend aus drei Stationen mit jeweils zwei EVU-Stationea Rg. 2 zeigt der , Aufbau 
Fe Wes! Reiner EVU-Station. der sich durch die vetwendeten Symbols und die Beschrrftung selbst erWart. 
Die ne ^Spologfeche Adresse ergibt sich dann als Datenstruktur. die fQr jede dieser H.erarch.estoifen emen 
IdenSLTenthTzXeeispie. kann I topologische Adresse des einen Endes der Leitung 1 .n Staton 1 (s.ehe 
Rg. 1 und 2) folgenderma&en spezifiziert werden: 
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Station: 


1 


2 


EVU-Station: 


1 


3 


Feld: 


10 


4 


Terminal: 


10 



und folgendermaBen gelesen werden: Die beschriebene topologische Adresse befindet sich in Station 1 , innerhalb der 
Station 1 in der EVU-Station 1 , innerhalb der EVU-Station im Feld 1 0 und innerhalb des Feldes 1 0 am Terminal 1 0. Fig, 
1 und Fig. 2 verdeutlichen die Vorgehensweise. 

Fur die Grenze jeder Hierarchiestufe gilt, daft diese nicht uberlappen durfen. Das heiftt, jedes Terminal darf nur zu 
einem Feld, jedes Feld nur zu einer EVU-Station und jede EVU-Station nur zu einer Station gehOren. 

Aus der technologischen Beschreibung der Netztopologie kOnnen nun die abstrakteh Endpunkte, d.h. die internen 
topologischen Adressen, des Netzgraphen Qber die folgend beschriebene Datenstruktur, die von der Datenaufbereitung 
bereitgestellt wird, sehr einfach und schnell ermittelt werden. 

Die Transformation der externen in die interne topologische Adresse ist in Rg. 3 dargestellt und wird nachstehend 
beispielhaft beschrieben: 

Zwischen EvuStalnStati6n[i} und EvuStalnStation[i+1] liegen alle EVU-Stationen der Station i. In der Datenstruktur 
EvuStalnStationO s\nd die Identifizierer enthalten, so daB die jeweiligen Teiladressen nicht dicht liegen mussen, sondern 
durch direkten Vergleich ermittelt werden. Ausnahme ist der Stationsidentifizierer, der als Index fur EvuStalnStationO 
dient. Entsprechend werden die anderen Hierarchiestufen durchlaufen. Der schlieBlich ermittelte Index der Datenstruktur 
TerminafQxsX die interne topologische Adresse des mit der externen topologischen Adresse spezif izierten Terminalpunk- 
tes. Die interne topologische Adresse wird als Index fur die Endpunkte (irrterner Terminalindex) des Netzgraphen ver- 
wendet. 

Durch die BerOcksichtigung der Felder als Teiladresse ist die Ermittlung der Feldtopologie als impliziter Teil der 
Netztopologie mdglich. Hierdurch kfinnen die Berechnungsfunktionen zur Ermittlung der lokalen Feldtopologien entfel- 
len. 

Auch MeBwerte und Schutzmeldungen lassen sich als topologische Objekte modellieren. Dies vereinfacht die 
Dateneingabe und ermOglicht die exakte Beschreibung aller Netzkonfigurationen, bei denen bisher zusatzlich Pseudo- 
Netzelemente vorzusehen waren. 

MeBwerte sind ferngemeldete analoge Werte relevanter ProzessgrdBen wie Spannung, Strom und Leistung etc., 
die im Leitrechner angezeigt und auf Grenzwertverletzung uberwacht werden. Bei Leitsystemen, die ebenfalls ein Topo- 
logiemodell im Leitrechner halten, werden die MeBwerte in Verbindung mit der Netztopologie und den elektrischen 
Parametern der Netzelemente zur Berechnung eines konsistenten und vollstandigen Netzzustandes verwendet (State 
Estimation). Hierzu mussen LeistungsmeBwerte jeweils den adjazenten impedanzbehafteten Netzelementen und Span- 
nungsmeBwerte den Knoten des Netzmodells zugeordnet werden. Als Knoten bezeichnet man hierbei alle impedanzlos 
zusammengeschalteten Netzelemente. Knoten werden durch impedanzbehaftete Zweige wie Lertungen und Transfbr- 
matoren verbunden. Bisher werden bei der Dateneingabe MeBwerte und Netzelemente statisch verknQpft. Dies hat den 
Nachteil, daBfQrbestimmte Netzkonfigurationen die MeBwerte andersplaziert werden mGssen, als es ihrertatsSchlichen 
Netzanordnung entspricht. Diesen gravierenden Nachteil vermeidet man, wenn LeistungsmeBwerte als gerichtete impe- 
danzlose Zweige und SpannungsmeBwerte als Objekte mit einer externen topologische Adresse modelliert werden. 
Die schaltungsabhangige Zuordnung von MeBwert und Netzelement wird dann von der Netztopologie geleistet. Wei- 
terhin wird die Zahlpfeilrichtung der MeBwerte direkt von der Orientierung des MeBwertzweiges festgelegt. In Rg. 4 ist 
die beschriebene topologische Zuordnung dargestellt. Mit P sind Wirkleistungswerte, mit Q Blindleistungswerte bezeich- 
net. Je nach Schaltzustand der Schafter S1 bis S4 ergeben sich die in folgender Tabelle aufgelisteten MeBwertzuord- 
nungen und deren Eigenschaften, wobei L1 , L2 Leitungen sind. 
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Tabelle 1: Topologische MeGwertzuordnung 



Wahrend bei der statischen Zuordnung der MeBwert nur einem Netzelement zugeordnet werden kann, ist durch 
die topologische Zuortnung in alien Fallen eine richtige Zuordnung bei gleichzeitig vereinfachter Dateneingabe mOglich. 
SinngemSB gilt die gleiche Vorgehensweise fur gerichtete Schutzmeldungen wie ErdschluBwischermeldungen und Kurz- 
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schluBanzeiger, die fur die Lokalisierung von ErdschlQssen oder Kurzschlussen von der Topologiefunktion ausgewertet 
werden. 

Unter Anwendung der vorstehenden Adressen-Definition erfolgt nachstehend eine ausfuhrliche beispielhafte 
Beschreibung des erf indungsgemaBen Verfahrens. Alie Verfahrensschritte lassen sich als Rechnerprogramm realisieren 
5 und mit Hilfe einer Datenverarbeitungseinrichtung automatisiert durchfuhren. 

1 . Auf bau der statischen Topologie 

Der Graph des betrachteten Netzes besteht aus Endpunkten, sogenannten Terminals, sowie impedanzlosen und 

10 impedanzbehafteten Zweigen (Zweitore) wie Schalter, MeBwerten, Verbindern sowie Leitungen, Transformatoren und 
LSngsblindleistungselementen. Weiterhin gibt es zusatzlich Eintore wie Generatoren und Lasten, SpannungsmeBwerte 
und Erdungsschalter etc. , die nur eine topologische Adresse sowie verschiedene Merkmale und Eigenschaften besitzen, 
die sie auf den zugehCrigen Teilgraphen Qbertragen, wie z.B. geerdet, spannungslos oder unter Spannung. 

Die Datenstruktur des Verfahrens basiert auf einer Terminal-Adjazenz-Matrix, die als vorwarts verkettete Liste abge- 

15 legt wird und bei der fur jedes Terminal die adjazenten Zweige zugeordnet werden (im folgenden mit Csftfbezeichnet). 
Fur diese symmetrische Matrix, die reine Strukturinfbrmationen enthait, wird eine optimale Eliminationsreihenfolge 
bestimmt und CstQ urn die hierbei zusatzlich entstehenden Fill In (f iktive Verbindungszweige) erganzt. Aufgrund der 
Ahnlichkeit mit der numerischen LDU-Zerlegung einer Matrix wird im folgenden dieser Schritt strukturelle Dreiecks- oder 
LDU-Zerlegung genannt. Unter LDU-Zerlegung wird die Zerlegung einer Matrix in das Produkt einer Lower-, Diagonal- 

20 und einer Upper-Matrix verstanden, wobei die Lower- und Upper-Matrix jeweils Dreiecksmatrizen sind. Urn die topolo- 
gische Eigenschaft der Fill In zu verdeutlichen, werden diese im folgenden als f iktive Zweige bezeichnet Wie im folgen- 
den gezeigt wird, enthait die strukturelle Upper-Matrix (im folgenden mit ParkettO bezeichnet) Information uber die 
Konnektivitat des Graphen, die bei Netzanderungen sehr einfach und schnell in ParkettO aktualisiert werden kann. 
Anhand des in Fig. 5 gezeigten Weinen Beispielgraphen werden die skizzierten Datenstrukturen verdeutlicht. Alle 

25 im Graphen abgebiideten Zweige sind Bestandteil der statischen Topologie. Die mit GroBbuchstaben gekennzeichneten 
Zweige haben den Status EIN, die mit Kleinbuchstaben gekennzeichneten den Status AUS. In diesem Beispiel sind nur 
die internen Terminalindizes und Zweigbezeichner aufgefuhrt. Die AbbikJung von externer auf interne Topologieadresse 
geschieht mit Hilfe der eriauterten Adresshierarchie und kann aus Ubersichtsgrunden hier einfach weggelassen werden. 
Ebenfalls ist der Zweigtyp in diesem Beispiel von untergeordneter Bedeutung. Auf eine detaillierte Beschreibung wird 

30 deshalb, ohne die Allgemeingultigkeit des Beispiels einzuschranken, verzichtet Das Beispielnetz mit 9 internen topo- 
logischen Adressen und 1 6 Zweitoren oder Zweigen wird in den folgenden Stuf en auf die oben skizzierte Datenstrukturen 
abgebildet. 

Die Terminal-Adjazenz-Matrix fOr diesen Graphen und deren LDU-Zerlegung kann nun wie in Fig. 6 notiert werden. 
Die Matrixelemente enthalten zur naheren Eriauterung des Algorithmus die Bezeichner der zugehdrigen Zweige. Die 

35 mit a gekennzeichneten Elemente der Matrix entsprechen den Fill In, d.h. den f iktiven Zweigen im Netzwerk. Die Rei- 
henfolge der Terminals ist bei diesem Beispiel zur Vereinfachung schon optimal gewahlt. Die interne Programmablage 
der Terminal-Adjazenz-Matrix geschieht als gedrangt abgespeicherte, vorwartsverkettete Liste (Csffl), die der LDU- 
Matrix als eindimensionales Array {ParkettO) von C-Strukturen. Wegen der Struktursymmetrie wird nur cfie Upper-Matrix 
abgelegt. Im folgenden wird angenommen das alle Zeilen der Upper-Matrix und alle Elemente dieser Zeilen in optimaler 

40 Reihenfolge sortiert sind. Jedes Objekt der verketteten Liste bzw. des Arrays enthait fur jedes Nichtdiagonal-Element 
den Terminalindex des gegenQberliegenden Terminals, den Index in der Zweigliste, den Status des Zweiges und den 
Zweigtyp. 

2. Initialisierung der Upper-Matrix (vergl. Fig. 7, 8 und 10) 

45 

Jedes Element der Upper-Matrix wird aufgrund des zugehdrigen Zweigstatus initialisiert: 
Status EIN => z.B. K = SCHAR_MIN (-128, minimaler Wert eines Byte mit Vorzeichen); Status AUS: 0. Die fiktiven 
Zweige erhalten den Status AUS (0) . FOr die so initial isierte Matrix wird folgende Prozedur durchgef Qhrt und in Anlehnung 
an die rein ganzzahligen Operationen ganzzahlige Initialisierung genannt. Diese geschieht analog dem Matrix- Produkt 
so aller Zeilen ab dem ersten Nichtdiagonalelement rechts der Diagonalen mit sich selbst: 



55 
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Fiir alle Zeilen i 

Fur alle Zeilenelemente j, k der Zeile i rechts der 
Diagonalen, wobei i=l..n, j=i+l..n, k=j..n und n^Anzahl 
Terminals: 

Falls die zugehorigen Matrixelemente und (k,i) 

ungleich 0 sind: 

inkrementiere die Elemente aller zugehorigen 
Schnittpunkte des Zeilen- und Spaltenbruchstuckes 
rechts bzw. unterhalb der Diagonalen in der Upper- 
Matrix urn 1. Durch die Struktursymmetrie 1st das 
Zeilen-und Spaltenbruchstiick identisch 

Fur die erste Zeile in Fig. 7 ist der erste nichtdiagonale 
Schnittpunkt das Matrixelement (5, 8) und diesem Element is 
der Zweig f zugeordnet. Dieses Element ist mit 0 
initialisiert und erhSlt den aktuellen Wert von 1. 



In Fig 7 ist dies for die erste Zeile veideutlicht-Nachdieserganzzahligen Initialisierung liegtnundiein Fig.8dargestellte 
Upper-Matrix vor. Normalgedruckte Kleinbuchstaben entsprechen ausgeschalteten Zweigen und haben den StatuszSh- 
ler 0 normalgedruckte GroBbuchstaben entsprechen eingeschaKeten Zweigen und haben den Statuszahler 
SCHAR MIN + i wobei i durch die Anzahl der auf dieses Element ausgefOhrten Inkremerrtoperationen bestimmt ist. 
Die mit a gekennzeichneten Elemente der Matrix entsprechen den Fill In, d.h. den fiktiven Zweigen im Netzwerk. die 
entsprechend der aktuellen Topologie aufgrund der ganzzahligen Operationen den Zustand EIN oder AUS erhalten. 
Dies entspricht einem Statuszahler 0 Oder i, wobei i durch die Anzahl der auf dieses Element ausgefuhrten Inkrement- 
operationen bestimmt ist. Die mit felt gedruckten Kleinbuchstaben bezeichneten Zweige haben zwar den aktuellen Sta- 
tus AUS, erhalten aber aufgrund der ganzzahligen Operationen den Status eines verbindenden Zweiges und haben 
ebenfalls einen Statuszahler von i, mit i = Anzahl Inkrementoperationen. Weiterhin kann der Status des Zweiges auch 
von seinem Typ abhangen, z.B. urn spezielle Topologien, wie z.B die Netzbezirkstopologie zu ermitteln. Zur ErmittJung 
der Netzbezirkstopologie sind alle Transformatorenalstrennende Zweige vorzugeben. 

Die Konnektivitat des Graphen laBt sich nun sehr einfach aus der Pfad-Tabelle (PathTableB) der Upper-Matrix Par- 

/cerf/7ableiten(vergl. Rg.9). . 

Diese Pfad-Tabelle errthalt fur jeden Terminalindex in ParkettU den gegenOberhegenden Terminalindex des ersten 
adiazenten Zweigs Diese Tabelle ist neben der Upper-Matrix die zweite wesentliche Datenstruktur dieses Topologie- 
verfahrens Aus ihrkanndie Konnektivitat des Netzgraphen leicht ermrttelt werden. Diese Tabelle enthalt so viele Emtrage 
mit dem Wert 0. wie der Netzgraph in Teilgraphen zerfallt Die in Fig. 9 dargestellte Pfad-Tabelle kann aus obiger Struktur 
Parte/fl abgeleitet werden. Diese Pfad-Tabelle enthalt nur eine 0 und der Netzgraph ist deshalb zusammenhangend, 
wie leicht am Pfad-Graph in Fig. 9 zu erkennen ist 

Durch Zuordnung z.B eines Farbbezugpunktes f Or die Einf arbung des Netzgraphen oder eines Merkmals zur Kenn- 
zeichnung bestimmter Zustande zu einem Entry dieser Pfad-Tabelle (i.e. interne topologische Adresse oder Terminal- 
index), kann nun dem entsprechenden Teilgraphen die zugehOrige Faroe oder Eigenschaft zugeordnet werden. indem 
der P/ad-Gr^p/7durchlaufen wird. 
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3. Update der Netztopologie(vergl. Rg. 10, 1 1 und 12) 

Angenommen im Beispielgraph dndert sich der Status des Zweiges G von E/A/nach AUS. Der ganzahlige Updateal- 
gorithmus ermittelt nun Anderungen for mit => gekennzeichnete Matrixelemerrte, wie Rg. 11 zeigt. 
5 Der Algorithmic kann fur eine Ausschaltung bzw. Einschaltung eines Zweiges folgendermaBen formuliert werden: 



Ausschaltung eines Zweiges : 

Subtrahiere vom entsprechenden Element der Upper Matrix die 
Konstante K. Falls hierdurch dieses Element zu 0 wird 1 : 

Update AUS: 

Fur dieses Element: 

dekrementiere fur alle Schnittpunkte dieses Elementes mit 
alien anderen Elementen ungleich 0 der zugehorigen Zeile 
bzw. Spalte rechts bzw. unterhalb der Diagonalen, die 
entsprechenden Elemente in der Upper-Matrix um 1. Fur 
jedes Element , das hierdurch zu 0 wird, wiederhole 
Update AUS 

Einschaltung eines Zweiges : 

Addiere zum entsprechenden Element der Upper Matrix die 
Konstante K. Falls hierdurch dieses Element zu K wird 1 : 

Update Ein : 

Fur dieses Element: 

dekrementiere fur alle Schnittpunkte dieses Elementes mit 
alien anderen Elementen ungleich 0 der zugehorigen Zeile 
bzw. Spalte rechts bzw, unterhalb der Diagonalen, die 
entsprechenden Elemente in der Upper-Matrix um 1 . Fur 
jedes Element f das hierdurch zu 1 wird, wiederhole 
Update EIN 

l Nur dann hat die Zustahdsanderung einen topologischen EinfluB 
und die Aktualisierung der Upper-Matrix darf stattfinden 



Fur dieses Beispie! fuhrt die Dekrementierung des Matrixelementes (8,9) dazu, da3 dieses Element den Wert 0 erhait. 
Daraus folgt, daft die aktualisierte Pfad-Tabelle einen weiteren Null-Eintrag erhait, d. h. der Netzgraph zeriaift in zwei 
Teilgraphen. Die in Rg. 12 dargestellte Pfad-Tabelle I33t sich aus der aktualisierten Struktur ParkettQ ermitteln. Es ist 
an der Pfad-Tabelle und am Pfad-Graph offensichtlich, daG der Graph in zwei Teilgraphen zerfallen ist. Eine Farbgebung 
55 oder die Weitergabe eines Merkmals ist durch die Bearbeitung der Pfad-Tabelle nun sehr einfach und sehr schnell 
mGglich. 
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4. LastfluBschrilttstelle 

Die Struktur der LDU-Matrix fur LastfluB- Oder KurzschluBberechnungen erhalt man in zwei Stufen. In der ersten 
Stufe werden alle impedanzlosen Zweige. auBer den MeBwerten, zu Knoten zusammengefaBt. AnschlieBend werden 
5 die MeBwerte den impedanzbehafteten Netzelemerrten zugeordnet und in der zweiten Stufe die MeBwerte ebenfalls 
mit den Knoten verschmolzen. Rg. 13 zeigt ein Knoten-Zweig-Netz. Unterstellt man, daB die Zweige d.e.f,G.K impe- 
danzlos sind. so kann aus ParkettQ unmitlelbar die (LD)U-Form der Knotenadmittanzmatrix extrahiert werden. Auf eine 
MeBwertzuordnung wird in diesem Beispiel verzichtet. 

Die Terminals 1 ,5,6,8 verschmelzen in den Knoten 8. wahrend alle anderen Terminals durch impedanzbehaftete 
10 Zweige verbunden sind und deshalb erhalten Weiben. Es ergibt sich also die in Rg. 14 dargestellte Upper-Form der 
strukturellen JACOBI- oder Knoten-Admittanz-Matrix. _ . , . 

Die Filter-Funktion elirriniert alle impedanzlosen Zweige und transferiert die betroffenen Zweigenden (Terminals) 
von impedanzbehafteten Zweigen zu den Terminals mit hOherem Index, wie hier z.B. 1.5.6-»8. Dies ist somit schon die 
Struktur der JACOBI-Matrix (Lastf luB) oder der Knoten-Admittanz-Matrix (KurzschluB). Die weiteren Schritte einer Last- 
is fluBrechnung sind dann einfach durchtOhrbar. Indem die LDU-Zerlegung schon auf der Ebene der Basis-Topologie 
(Schalter-Terminal) durchgefuhrt wird. kOnnen eine Vielzahl von Programmschritten entfallen. Diese Vereinfachung der 
Programme verbessert im erheblichem MaBe die Rechenzeit und die Fehleranfalligkeit der Berechnungsprogramme 
wie State-Estimation, LastfluB und KurzschluB. 

20 5. Testnetze 

Rg. 15 zeigt die Form eines generierbaren Netzes, mit dem Performancetests for den Algorithmus durchgefuhrt 
wurden. Durch Parametervorgabe kann der Vermaschungsgrad VGdes Netzes zwischen 0 und 3 und die NetzgrOBe 
beliebig vorgegeben werden. Der Vermaschungsgrad ergibt sich aus (1). 
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VC- Anzahl Zwe '9 e (1) 
Anzahl Terminals-! 



In elektrischen Obertragungsnetzen und Verteilnetzen variiert der Vermaschungsgrad zwischen 1 und 1 ,5. Fur ein gro- 
Bes reales Verteilnetz ermittelt sich z.B. der zugehorige Vermaschungsgrad VQ inWusive der einspeisenden 11 OkV 
30 Stationed und Umspanner folgendermaBen: 





Anzahl Schalter (ohne Efdungsschalter): 


9683 


35 


Anzahl Leitungen: 


2100 




Anzahl 110/10kV-Umspanher: 


22 




Anzahl 10/0,4KV-Nelztrafos: 


1200 


40 


Anzahl Zweige: 


13005 



Mit der Anzahl Terminals von 1 1 .623 ergibt sich ein Vermaschungsgrad von VG = 1 ,11 9. Das bedeutet. daB die im 
folgenden vorgelegten Performancezahlen for VG = 1 .2 bei groBen Verteilnetzen als Obergrenze anzusehen sind. 

6. Performance 

6.1. Anzahl fiktiver Zweige Qber die NetzgrdBe in AbhSngigkeit vom Vermaschungsgrad 

Wie oben beschrieben, wird die Terminal-Adjazenz-Matrix in eine obere Dreiecksmatrix transformiert. Hierbei war- 
den fiktive Verbindungszweige (Fill In) erzeugt. Die Anzahl zusatzltoher Verbindungszweige ist abhSngig von der Rei- 
henfokje der einzelnen Transformationsschritte und vom Vermaschungsgrad VG. Beim Aufbau der oberen 
Dreiecksmatrix wird die Reihenfolge der Transformationsschritte optimiert. Deshalb reduziert sich die Abhangigkeit auf 
den Vermaschungsgrad und die Anzahl Zweige. wie im Rg. 16 dargestellt. Die Anzahl fiktiver Zweige beeinfluBt die 
Performance des vorgeschlagenen Verfahrens. Die weiteren Diagramme fOr die einzelnen Programmschritte zeigen 
dies deutlich. 
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6.2. Performance der einzelnen Programmschritte Liber die NetzgroBe in AbhSngigkeit vom Vermaschungsgrad 

Die CPU-Zeit eines 80486 40MHz Rechners wurde fur folgende Programmschritte in Abhangigkeit der Anzahl 
Zweige und des Vermaschungsgrades gemessen 

1 . Aufbau der slatischen Topologiestruktur ParkettO 

2. Initialisierung der dynamischen Topologie auf Basis der aktuellen Status aller Netzzweige 

3. Update der dynamischen Topologie aufgrund einer Statusanderung eines Zweiges 

Fig. 1 7 zeigt die benOtigte CPU-Zeit in 1 0 Mips*s fur den Aufbau der statischen Topologie ParkettO- Die Zeit fur eine 
100 Mips CPU ergibt sich dann durch Division mit 100 Mips. 

Der Aufbau dieser Struktur ParkettQ\gt der aufwendigste der 3 Verarbeitungsschritte. Ein Aufbau dieser Struktur ist 
allerdings nur bei Netzausbau notwendig. Legt man eine 200 Mips CPU (z.B. DEC AXP 600) zugrunde, so benfitigt man 
f Or ein Netz ahnlicher GrOBe wie das unter Punkt 5 beschriebene Testnetz ca. 3s CPU-Zeit. 

Die GrOBe der Struktur Parfreff/7 wird folgendermaBen berechnet: 
(2*Anzahl Zweige+Anzahl Terminals)*16 Byte 
Fur ein Netz in der GrOBenordnung des Testnetzes ergibt sich die GrOBe von ParkettO zu ca. 0.6 MByte. 

Fig. 18 zeigt die bendtigte CPU-Zeit in 10 Mips*s fur die Initialisierung der statischen Topologie ParkettO mit den 
aktuellen Status aller Netzzweige: 

Eine Initialisierung der Topologie ist typischerweise beim Systemstart, nach Datendnderungen und je nach Doppelrech- 
ner-Philosophie beim Rechnerwechsel notwendig. Fur ein Netz in der GrOBenordnung des Testnetzes ware die Initiali- 
sierung auf einer 200 Mips CPU in ca. 25 Millisekunden CPU-Zeit abgeschlossen. 

Fig, 1 9 zeigt die benOtigte CPU-Zeit in 10 Mips*ms fur den Update der Topologie ParkettO mit dem aktuellen Status 
eines Netzzweigs: 

Die mittlere Zeit fur einen Update der Topologiestruktur bei einer Statusanderung z.B. eines Schalters betragt fur ein 
Netz in der GrfiBenordnung des Testnetzes fQr eine 80486 CPU mit 40MHz ca. 8.0 10"5 s. FOr eine 200 Mips CPU liegt 
sie in der GrOBenordnung von ca. 4.0 10~ 6 s. Diese Zeit ist urn mehrere GrOBenordnungen schneller, als die Zeiten der 
momentan eingesetzten und bekannten Verfahren. Hierdurch sind ganzlich neue Konzepte bei der Topologieverarbei- 
tung sowie bei den anschlieBenden Berechnungsverfahren mOglich. Zum Beispiel kOnnen alle Zweige des Netz es nach 
jedem Update darauf untersucht werden, ob der jeweilige Zweig die einzige Verbindung (kritischer Zweig) zwischen 
zwei Teilnetzen ist. FQr ein Netz in der GrOBenordnung des Testnetzes und eine 200 Mips CPU wurde die CPU-Gesarrrt- 
zeit inclusive Update ca. 

2*Anzahl Zweige *T updatQ * 100 Millisekunden 

dauern. Weiterhin zeigt dieses Beispiel, daB die kumulierte CPU-Zeit fur die Ausfuhrung von Updates fQr alle Scharter 
in diesem Beispiel ca. 50 Millisekunden betragt Dies ist nur etwa die doppelte Zeit, die die Initialisierung bendtigt und 
ein weiterer Hinweis, daB der Update bei diesem Verfahren nur die wirklich notwendigen Anderungen der Datenstruktur 
umfaSt. 

Patentanspruche 

1. Verfahren zur Initialisierung und Aktualisierung eines Netzmodells, wobei die Netztopologie in Form eines Netzgra- 
phen unter Anwendung und Erweiterung von Sparse- Vector-Verfahren und Sparse-Matrix-Verfahren durch nach- 
stehende Schritte gebildet wird: 

a) Bereitstellung erfaSter Netztopologiedaten in Form einer Terminal-Adjazenz-Matrix; 

b) Umformung der Terminal-Adjazenz-Matrix in eine obere Dreiecksmatrix, wobei zusatzliche fiktive Verbin- 
dungszweige entstehen kOnnen; 

c) Initialisierung der oberen Dreiecksmatrix, wobei jedes Element der oberen Dreiecksmatrix aufgrund des zuge- 
hOrigen Zweigstatus initialisiert wird, und wobei 

c1) jede zu einem eingeschalteten Zweig gehdrendes Element mit einer Weinen negativen ganzzahligen 
Konstante K initialisiert wird, die so gewahlt ist, daB sie durch die weiteren Verarbeitungsschritte nicht 
grOBer als Null werden kann, und 

c2) jede zu einem ausgeschalteten oder fiktiven Zweig gehdrendes Element mit Null initialisiert wird, 

d) Bildung der Matrixprodukte fur alle Zeilen i (i = 1 bis maximale Anzahl der Zeilen), der oberen Dreiecksmatrix 
ab und ausschlieBlich der Diagonalen nach folgender Regel: falls die Elemente i,j und i,k ungleich Null sind, 
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sind die Elemente atler Schnittpunkte des Zeilen- und Spalten-BruchstOckes der oberen Dreiecksmatrix urn 1 
zu inkrementieren, wobei j und k Laufindizes fur alle Elemente der betreffenden Zeile sind; 

e) Bildung einer Pfad-Tabelle aus der oberen Dreiecksmatrix. wobei die Pfad-Tabelle fQr jede Zeile bzw. jedes 
Terminal der oberen Dreiecksmatrix den inzidenten Terminal-Index des ersten von Null verschiedenen Elemen- 
tes dieser Zeile angibt; 

f) Aktualisierung der oberen Dreiecksmatrix im Fall einer Anderung der erfaBten Netztopologiedaten, also einer 
Statusanderung eines Netzzweiges, durch folgende MaBnahmen: 

f1) im Fall einer Ausschaltung eines Netzzweiges: vom entsprechenden aktuellen (= bisherigen) Wert des 
Matrixelementes wird die Konstante K subtrahiert; 

f1 1) falls dieses Element dadurch zu Null wird, sind die Elemente aller Schnittpunkte der oberen Drei- 
ecksmatrix dieses Zeilen- und Spalten-Elementes mit alien anderen Elementen ungleich Null derzuge- 
hOrigen Zeile bzw. Spalte. rechts bzw. unterhalb der Diagonalen urn 1 zu dekrementieren; fails dadurch 
weitere Elemente zu Null werden, wird fQr diese die Operation in f 1 1) wiederholt; 

f2) im Fall des Einschaltens eines Netzzweiges: zum aktuellen Wert des Matrixelementes wird die Konstante 
Kaddiert; 

f21) falls dieses Element dadurch zu K oder fur weitere Wiederholungen zu 1 wird, sind die Elemente 
aller Schnittpunkte der oberen Dreiecksmatrix dieses Zeilen- und Spalten-Elementes mit alien anderen 
Elementen ungleich Null der zugehOrigen Zeile bzw. Spalte, rechts bzw. unterhalb der Diagonalen urn 
1 zu inkrementieren; falls dadurch Elemente zu Eins werden, wird fQr diese die Operation in f21) wie- 
derholt; 

g) Aktualisierung der Pfad-Tabelle entsprechend dem Schritt e); 

h) Kennzeichnung der Komponenten des Netzgraphen anhand der Pfad-Tabelle, wobei ausgehend vom ersten 
Terminal-Index der Tabelle alien Terminals des zugehCrigen Teil-Pfades die gleiche Komponenten-Kennzeich- 
nung zugeordnet wird; ergibt sich beim Durchlaufen weiterer Teil-Pfade eine schon vergebene Kennzeichnung 
(durch Treffen auf eine Kbmponente, der schon eine Kennzeichnung zugeordnet ist), so wird diese Kennzeich- 
nung dem gerade durchlaufenen Teil-Pfad zugeordnet; der Netzgraph zerfailt in Teil-Graphen entsprechend der 
Anzahl vergebener Kennzeichnungen. 

Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB die Konstante K eine ausreichend groBe positive Zahl 
ist, und anstatt Dekrement- Inkrementoperationen und umgekehrt auf die Matrixelemente ausgefQhrt werden. 

Verfahren nach Anspruch 1 Oder 2, dadurch gekennzeichnet, daB zur Bereitstellung der Netztopologiedaten in einer 
zweckma&igen Datenstruktur interne und externe topologische Adressen verwendet werden, wobei ausgehend von 
einer definiertentopologischen Hierarchie jede Hierarchiestufe miteinem Identifiziererbezeichnetwird, und in einer 
externen topologischen Adresse eines Terminals (niedrigste Hierarchiestufe) die ZugehOrigkeit zu bestimmten Kom- 
ponenten der hOheren Hierarchiestufen enthalten ist, woraus die abstrakten Endpunkte des Netzgraphen als interne 
topologische Adressen ableitbar sind. 

Verfahren nach Anspruch 3, dadurch gekennzeichnet, daB MeBwerte und Schutzmeldungen topologisch zugeord- 
net werden, wobei im Fall des Modells eines elektrischen Netzes LeistungsmeBwerte als gerichtete impedanzlose 
Zweige und SpannungsmeBwerte als Objekte mit einer externen topologischen Adresse modelliert werden. 
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